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Abstract. The Cauchy problem for the Vlasov-Maxwell-Boltzmann equations (VMB) 
is considered. First the renormahzed solution to the Vlasov equation with the Lorentz 
force is discussed and the difficulty on the partial differentiability of the coefficients is 
overcome. Then the weak stability of the renormalized solutions to the Cauchy problem 
of VMB is established using the compactness of velocity averages and a renormalized 
formulation. Furthermore, the large time behavior of the renormalized solutions to VMB 
is studied and it is proved that the density of particles tends to a local Maxwellian as 
the time goes to infinity. 



1. Introduction 

Since the work of DiPerna and Lions [lOj on the Cauchy problem for the Boltzmann 
equation twenty years ago, it has been a weh-known open problem to extend their the- 
ory to the Vlasov-Maxwell-Boltzmann equations. Among the difficulties, how to define 
the characteristics of the Vlasov-Maxwell-Boltzmann equations is a major obstacle. In 
this paper, we will give the following partial results: the weak stability and large time 
behavior of the renormalized solutions to the Vlasov-Maxwell-Boltzmann equations, and 
existence of the renormalized solutions to the Vlasov equation with the Lorentz force. The 
fundamental model for dynamics of dilute charged particles is described by the Vlasov- 
Maxwell-Boltzmann equations (VMB) of the following form [51 [3 HSl [HI [231 127] : 

^ + ^-^xf + (.E + (xB)-V^f = Q{f,f), xeW", CeM^, t>0, (1.1a) 

_ - V X 5 = -J, divB = 0, on x (0, oo), (1.1b) 
dB 

— + VxE = 0, divE = p, on Mi^x(0,oo), (1.1c) 

p= [ fdC, j= [ fCd^, on M3x(0,oo), (l.ld) 

where / = /(t,x,^) is a nonnegative function for the density of particles which at time t 
and position x move with velocity ^ under the Lorentz force 

E + ^xB, 
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E is the electric field, B is the magnetic field, the function j is called the current density, 
and the function p is the charge density. The collison operator Q{f,f), which acts only 
on the velocity dependence of / (this reflects the physical assumption that collisions are 
localized in space and time), is deflned as 

Q{fJ)= [ dojb{C-C*,uj)if'f:-ff*), 

Jr3 Js2 

with (jj E 5^, the unit sphere in M^, where b = b{z,uj) denotes the collision kernel which is 
a given nonnegative function deflned on x 5^, and 

f* = f{t,x,C.), f' = f{t,x,0, fi = f{t,x,C), 

with 

e' = e-(e-e*,^v, 

C = C* + {C- 

which yield one convenient parametrization of the set of solutions to the law of elastic 
collisions 

\e\' + \a' = \c\' + \^*\'- 

The interpretation of ^, (,1 is the following: are the velocities of two colliding 

molecules immediately before collision, while are the velocities immediately after the 
collision. Those unknown functions /, E, and B are strongly coupled, and the constraint 
on the divergence of E will be ensured provided that the conservation of charge holds; 
that is, 

dp , 

— + divr^j = 0, 

since 

d 

= —{diVxE - p) = diY^Et - pt 
at 

= div^(Va; X B - j) - Pt 
= -pt - div^j, 

due to the fact div(V x v) = for any vector- valued function v. Similarly, the magnetic 
fleld B remains divergence free if it is so initially. 

The VMB equations are integro-differential equations which provide a mathematical 
model for the statistical evolution of dilute charged particles. The construction of global 
solutions to VMB has been open for a long time until only a few years ago. In Guo |18) . 
a unique global in time classical solution near a global Maxwellian (independent of space 
and time) was constructed. See also Strain [27] for the extension to the Cauchy problem. 
Notice that. Lions constructed in [23] a very weak solution to VMB, which is usually called 
a measure-valued solution, using Young's measure to deal with the nonlinearity. 

For the particles without collision (cf. [H |9l [I5l [161 EH |26j ) , or when the molecules are so 
rare that they do not interact with each other, VMB becomes the so-called Vlasov-Maxwell 
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system (VM), 

^ + e-V,/ + (S + ex5)-V^/ = 0, x£R^ ^eW", t>0, (1.2a) 
OE 

_ - V X 5 = - j, divB = 0, on x (0, oo), (1.2b) 
dB 

— + \/xE = 0, diYE = p, on M^x(0,oo), (1.2c) 

P= I fdC, j= [ fCd^, on ]R3x(0,oo). (1.2d) 
Note that (jl.2ap is a transport equation with a divergence free coefficient, that is 

div,,^(e,^ + ex5) = o. 

This property ensures that the solution will remain the same integrability as the initial 
data. With the help of this observation and velocity averaging lemma, DiPerna and Lions 
proved in [9] the global existence in time of weak solutions to VM with large initial data. 
For the smooth solutions to VM, we refer the readers to Glassey |16) and Schaeffer |26j . 

The main goal of this paper is to show the weak stability and the large time behavior 
of the renormalized solutions to VMB. To this end, we will need an existence result of the 
renormalized solution to the Vlasov equation ()1.2ap . Notice that the Vlasov equation is 
a transport equation with only partially W^^^ regularity, since usually we can not expect 
any differentiability on the magnetic field B and the electric field E from the conservation 
of energy. Inspirited by the result in Bouchut [3J and Le Bris-Lions [20] , we will first show 
the existence of renormalized solutions to the Vlasov equation. The presence of a non- 
trivial magnetic field B(x,t), a natural consequence of the celebrated Maxwell theory for 
electromagnetism, creates severe mathematical difficulty in studying the weak stability of 
weak solutions and the construction of global in time solutions for VMB. Our first result 
on weak stability is built on our above mentioned new result about renormalized solutions 
to the Vlasov equation with the aid of the velocity average lemma (DiPerna-Lions [9] and 
DiPerna-Lions- Meyer [13] ) and some techniques from Lions |22^ [23] . Our second result on 
renormalized solutions to VMB is their large time behavior, since from the physical point 
of view, the density of particles is assumed to converge to an equilibrium represented by 
a Maxwellian function of the velocity as the time t becomes large. Our results heavily 
depend on, apart from the weak compactness property, 

• the existence of renormalized solutions to the Vlasov equation; 

• a renormalized formulation, which is crucial to make sure that the quadratic term 
Q{f, f) is meaningful in D' (sense of distributions); and 

• the velocity averaging lemma [U [13] , which is crucial for the convergence of nonnlin- 
ear term {E + ^ x B) ■ V^/. 

The stability of renormalized solutions under weak convergence yields a consequence 
on the propagation of smoothness for those solutions. Indeed, a sequence of renormalized 
solutions {fn}^=i to VMB is relatively strongly compact in L^([0,r] x M®) if and only if 
the sequence of the corresponding initial data {/onjj^i is relatively strongly compact in 
L^(M®). In other words, under our assumption on the collision kernel and the integrability 
of the electric field and the magnetic field, no oscillations develop unless they are present 
from the beginning. 
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In order to prove our results, the standard a priori estimates derived from the conser- 
vation laws and H theorem are very useful, and in addition we need some assumptions on 
the integrability of the electric field E{t^ x) and the magnetic field B{t, x). More precisely, 
besides the standard estimate of E and B in L°°(0, T; L^(M^)), we need to assume that E 
is uniformly bounded in L°°(0, T; ^^(M^)) and B is uniformly bounded in L°°(0, T; L^(R3)) 
for some s > 5. The reasons for these requirements on E and B are twofold: (I) 
when we define the characteristics for the Vlasov equation, we need a bound on E in 
L°°(0, T; L^(M^)); (II) the averaging lemma (cf. [23]), combining with the uniform bound 
of ^3 fdi in L°°(0, T; lI (M^)) and the uniform bound of ^3 ifd^ in L°°(0, T; lI (M^)), im- 
plies the compactness of the first two moments of / on U'{Q, T; L^^^{E?)) for any 1 < p < |, 
which is enough to ensure the convergence of the nonlinear Lorentz force term in the sense 
of the distributions provided that E and B are uniformly bounded in L°°(0, T; L^(]R^)) 
and L°°(0,r;L^(M3)) for some s > 5. 

We now remark that throughout this work we never claim the existence of renormalized 
solutions to VMB. Actually, all results in this paper are based on the assumption of such 
an exact existence or the existence of a sequence of approximating solutions. One possible 
direction to address the existence problem may be based on the construction of a sequence 
of exact solutions or approximating solutions with the requirement that the electric field 
E is uniformly bounded in L°°(0, T; L^(R^)). We notice that the hyperbolic property 
of the Maxwell equations also demonstrates some difficulties if we want to improve the 
integrability of the electric field and the magnetic field. How to fulfill this strategy is still 
an open question and will be the topic of our future research. 

When the Lorentz force disappears, that \s E + ^ x B = 0, VMB becomes the classical 
Boltzmann equation. For the Cauchy problem of the classical Boltzmann equation, in 
|10j DiPerna and Lions proved the global existence of renormalized solutions with angular 
cut-off collision kernel and arbitrary initial data, see also [U El El [9l |TT1 [22l [23] and the 
references cited therein. Later, Hamdache extended this existence result to a bounded 
domain in [19]. The method explored for the existence result was the analysis of the 
weak stability of solutions. The argument strongly relied on some compactness properties 
(see [22]) which hold for sequences of renormalized solutions. In [23], Lions extended 
the similar weak stability and global existence result to the Vlasov-Poission-Boltzmann 
equations. For the extension to the Landau equation, see Villani [28]. For the long time 
behavior of the Boltzmann equations, see [71 111 ! IT2 l I28j. 

This paper will proceed as follows. We will discuss the renormalized solution to the 
Vlasov equation in Section 2. Section 3 is devoted to stating a priori estimates for VMB, 
main assumptions and main results on the weak stability of renormalized solutions to 
VMB. Then, Theorem 13.11 on weak stability and Theorem 13.21 on the propagation of 
smoothness will be proved in Section 4 and Section 5, respectively. In Section 6, we 
study the large time behavior and establish mathematically the convergence of / to a lo- 
cal Maxwellian satisfying the Vlasov-Maxwell equations. Finally, in Section 7, we explain 
an extension of our results to the relativistic Vlasov-Maxwell-Boltzmann equations. 
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2. Renormalized Solutions to the Vlasov Equation 

In this section, we consider the Vlasov equation of the form: 

dtf + ^-VJ + {E + ^xB)-V^f = 0, (2.1) 

with B{x,t) G L°°(0,r;L2(M3)) and E{x,t) G L^{0,T; L"^ n L^{Rl)). 
If we set 

then (|2.ip becomes a standard transport equation 

dtf + B-Vf = 0. (2.2) 

The question of whether the Vlasov equation has renormalized solutions is not only useful 
when the normalized solution to VMB system is considered, but also has its own inter- 
est due to the lower regularity of the coefficients. The renormalized solutions mean that 
(j2.2p still holds if we replace / by f3{f) with a suitable /3. Over past twenty years, there 
are many important progress about the renormalized solutions to ()2.2p . More precisely, 
DiPerna and Lions showed in [5] the existence of renormalized solutions when the co- 
efficient B £ VF^'^(M^). In 2004, Ambrosio extended the DiPerna-Lions theory to BV 
(bounded variations) field in [2] (for related work, see [3]). Also, in 2004 Le Bris and Lions 
extended in [20] the DiPerna-Lions theory to the case that the coefficient has only partial 
regularity. 

For the VMB or the Vlasov equation, the velocity B is no longer in T^(^'^)^^^. In- 
spired by [21 [20] , we claim that we still can prove the existence of a renormalized solution 
to ([H]) under the conditions that E{x,t) G L°°{0,T; L^{R^) n ^^(M^)), and B{x,t) G 
L°°(0, T; L^(]R^)). This is a crucial step for establishing renormalized solutions to the 
Vlasov-Maxwell-Boltzmann equations. 

Theorem 2.1. Assume that B{x,t) G L°°(0, T; L2(IR3)) ^nd E{x,t) G L'^{0,T;L^n 
L^iRl)). Let fo G nL~(M6) and j^p/o e L^R^). Then there exists a solution to (l2T]) 
(and hence to (|2.2p ) such that 

f(t,x,0 e L^{[0,T],Ll^nL^^^{R^)), 

and I^P/ G L°°(0,T; L^{R^)), satisfying the initial condition f\t=o = foix,C)- Further- 
more, if fo G L^(L^(M3)), then f G L°°{0,T; L'^iL^iR^))), and hence the solution is 
unique. 

To begin with the proof, notice that B = {Bi,B2) satisfies 

diva;^! = div^;B2 = 0, 

with 

Bi{x,^) = ^ G W^f^^{R^) (it does not depend on x), 

B2ix, ^) = E + ^xBe LlioM^ <zL(^'))- 
The proof of this theorem is divided into three steps. The uniqueness is a crucial issue 
which is the consequence of the following two lemmas, the first one dealing with regular- 
ization, and the second one stating the uniqueness. Finally, we will show the existence 
part. 
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Now we denote the mollifier Kg clS 
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K G V{m?), / K = 1, K > 0, 

where = Co°(M'^). Then, we have the following two lemmas. 

Lemma 2.1. Let f = f{t,x,i) G L°°([0, T], Lf^^^^ nL^^^)(]R6)) he a solution of and 
Ks and K,^ be two regularizations with two different scalings, respectively, in the variable x 
and ^. Then, for any e > 0, there exists a number fx(e) with < fj,{e) < such that 

is a smooth (in (x,^)) solution of 
with 

lim^, = 0, ^n {[0,T], L^^^^^.^^n L^^^^^^^JR^)) . 

Lemma 2.2. Let f = f{t,x,S,) G L°° {[0,T], Lj^ ^-^ nL^^^(M^)) be a nonnegative solution 

of (USD with zero initial data fo = 0. //, in addition, e L'^{[0,T], L'^{M.^)) and 

f G Ll{L^), then / = for all time t > 0. 

We now prove these two lemmas, and then finally complete the proof of Theorem I2.1i 

2.1. Proof of Lemma 12.11 We will use the mollifier to regularize the function / in ^ 
and x, while we assume that / is differentiable with respect to t (the results below hold 
also for the general case from a standard mollification in t with the help of Lebesgue's 
dominated theorem.) All the functional spaces used here are local, which is clearly enough 
for such a regularization result. 

We first regularize in the variable by convoluting (j2.2p with to get 

+ (e • V,/) *K^ + {{E + CxB)- Vg/) * = 0. (2.3) 

Denoting by 

[{E + ^xB)- Vg, K^]{f) = {E + CxB)- V^if * K^) -K^*{{E + ixB)- V^/). 
Then, (12. 3p can be rewritten as 

+ (e • V,/) *n^ + {E + ixB)- Vg(/ * K^) = [{E + ixB)- Vg, K^]{f). (2.4) 

It is a standard fact (see [8J) that 

li^:=[{E + ixB)-V^,K^]{f)^Q in Lf,^^)^,,, (2.5) 

as — >■ 0. Indeed, this is clear for smooth coefficients and /, while the general case follows 
as in [8] by dense property through the estimate 

\\[{E + ixB)-V^,K^]{f)\\L. <C\\E + ixB\\^,,. WfU^, 

which then implies the following standard estimate by integrating in x, 

\\[{E + ixB)-V^,K,,\{f)\y <C\\E + ixB\\^, 1,1 J|/||l^- . 
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Next, we regularize in the x variable by convoluting (j2.4p with Kg for = / * to 
obtain, 



dfti * Ke 



dt 



(2.6) 



We now successively deal with each terms on the right-hand side of (j2.6p . First, it is 
easy to observe that for fixed ^, we have 

if * — > if, as e — ^ 0, 



in Lj^ j^^, which together with (j2.5p implies that 

lim lim * = 0, in LL 



(x',5),ioc" 



(2.7) 



/I— >-0 e— s>0 

Second, for the first term on the right-hand side of ()2.6p . we have 

[{E + ^y.B)- Vg, Ac,](^) = {E + iy.B)- V^{k, * f^) - * {{E + ^ x B) ■ V^/^) 

= {E + ^xB)- ((V^/^) * K,) - * ((i^ + e X B) • Vg/^) 
= [(ii; + exi?),K,](V^^). 

The latter bracket can be controlled as follows: 

[{E + ^XB),K,]{Vi:f^ 

Hence, for fixed /x, we have 

lim\{E + CxB)-V^,K^]if^) = 0, 

e->-0 

i'^ ^UoMc- '^^i^ implies, 

lim lim [{E + CxB)-V^, k^] (/^) = 0, (2.8) 

in L^^ j^^. By a standard diagonization procedure, for any e > 0, we can find fi{e) with 
< n{e) < e"^ such that 



<C\\E + ^XB\U W^MLr^.y 

(X ,4) ,LOC 



lim If * = 0, in LL 



£->-0 



(a;,^),/oc' 



and 



lim[(i? + Cxi?).V5,/t,.(,)](4(,))=0, in 



£^0 



To complete the proof of this lemma, it remains to show the following convergence for 
the above chosen fi{e): 



in L^^ Indeed, we can control /^^^-^'^ as 
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< 



|eV(;Ke|K^(e)|/(x - C,C " ri)\dCdr]. 



Thus, we deduce that, for any compact subset K C M^, by Fubini's theorem. 



< c 

< c 



K 



V ■ V(Ke(C)«^/,(e) {v)f{x - C, ^ - v)d(dr] 



e 



dxd^ 



eV^KeldC) sup \\f{x-C,C-r])\\L^ 
' ICI<£,l'7l<Me) 



(2.9) 



sup ll/(a;-C,C-??)llLi(jt)- 



^ ICI<e,l'7l<M{e) 

Since / G L|^^ ^-j, one has, according to the continuity of translation in L^[K), 

sup \\f{x - C,^ - - f{x,0\\L^K) ^ 0, as £^0, 

K\<e,\v\<Ke) 

and hence, 

sup — CtC ~ 'n)\\L'>^{K) is uniforly bounded for all e < 1. 

\(\<eM<K£) 

Thus, if we let e — ^ and < /u(e) < e^, we deduce from (j2.9p that 



I! 



/i{e),e 



0, in Ll 



as £ ^ 0. 



(2.10) 



Therefore, the lemma follows from (|2.7p . (12. Sh . and (I2.10p . and we complete the proof 
of this lemma. □ 



Next, we turn to the proof of Lemma [ 

2.2. Proof of Lemma 12.21 Let / be a nonnegative solution as claimed in Theorem 12.11 
We introduce two cut-off functions, respectively, with respect to each variable x and ^. 
For m,rz G N, we denote them by 



ipm{x) = ^ , and 0„(O = (t>(^ 



n 



where t/^ € P(M'^), tp > 0, tp = 1 for |x| < 1 and ip = for \x\ > 2; and the analogous 
properties are required on (j) with respect to the variable S^. We first multiply (j2.ip by (pn 
and integrate over ^ space to obtain 

4/ 



dt 



[ f(pndi + div, ( I U<Pndd\ + [ (E + ^xB)- V^f^ndC = 0. (2.11) 
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For the last term in (j2.1ip . we deduce, due to div^(£' + ^ x i?) = 0, 



/ (E + ^xB)- V^fckndi = - [ f{E + ixB)- Vi-cj^n 



Now we multiply (|2.1ip by V'm and integrate over x space to deduce 
d 



dt 



fiJm<t>ndxdS, + / IpmdiVr, ( / (,f(t>ndS, ) dx 

Jr3 VJrs / 



.l + \^\E + ^xB 



— 1 dS^dx. 
n 



(2.12) 



n 1 + lCl 

Hence, using the integration by parts for the second term in (j2.12p . we have 



d_ 
lit 



f'lprn(t>ndxd^ - / V^V' 



if(t>ndi dx 



(2.13) 



Next, we proceed to control the two integral terms in ()2.13p . Indeed, for the second 
term in (|2.13p . we have 



^xtpm ■ / if(t>ndi dx 



1 ^ , / ^ 
— Va;V' — 

3 m \m 



Cf(t>nd^ dx 



< 



1 



(2.14) 



C — \\^fX{m<\x\<2m,\i\<2n}\\L 



m 



as m ^ oo and n ^ oo. Here we used G L°°([0, T], (M^)), because 



\C\fdCdx < R 
< R 



fd^dx 



r3 7R3n{|5|>R} 



md^dx 



1 



fd^dx 



1 

R 



^\>R} 



\i?fd^dx 



<R\\f\\LHm + j^m?f\\L^ 



< 211/111. 



(R6)III^P/llll(R6) 



by optimizing the value of R. 

On the other hand, for m fixed, we claim that the term on the right-hand side of (I2.13P 
goes to zero as n goes to infinity by Lebesgue's dominated convergence theorem. Indeed, 
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as Vcf) is and supported in the annular {1 < |,^| < 2}, we have for almost all x E 



n 1 + |?| \n 



f 



n 1 + iei 



n 



< 2||V0||L-Vm||/X{n<|€|<2n}||Li(|^| + 1^1) 



^0, 

as n — )• oo, since for almost all x G 
Schwarz inequality, we have, 



fix,-) G 



In addition, by the Cauchy- 



f 



n 1 + 

l + \^\\E + ^xB 



f 



d^ 



n 1 + 1^1 ' 
<2\\VnL^\\fhi{\E\ + \B\) 

<i\m\L^i\\f\\l + \E\' + \B\')- 

and the right-hand side is in L],, since / E L^(L^) and E,B E L^. Thus, Lebesgue's 
theorem applies and we get the convergence of the term on the right-hand side of (I2.13P 
to zero as n goes to infinity, and m being kept fixed. 

Collecting the behaviors of those two terms, we obtain with (12.131) . as n, and next m, 
go to infinity, 

^ fdxdC = 0. 



□ 



dt 

As /o = 0, this yields / = for all t since / > and this concludes the proof. 



Having proved Lemma 12.11 and Lemma 12.21 ^ire now ready to complete the proof of 
Theorem 12.11 as follows. 

2.3. Proof of Theorem 12.11 Assume for the time being that we have at hand two 
solutions /i and /2 to (|2.ip satisfying the regularity stated in Theorem 12.11 and sharing 
the same initial value. In view of the interpolation between and L°°, and the fact 

E ( [0, T] , n L-^ (E3) ) n ( [0, T] , (M^ , (M^ )) ) 
for i = 1, 2, we deduce that 

heL°°[[0,T],Ll{R^LliR^))). 
By virtue of Lemma 12. H their difference f = /i — /2 satisfies 

^^^ + e-V(,,5)f^(,),, = A, (2.15) 
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with the same notation as in Lemma |2.1[ Since f /^(e) e ^ C°°(M^), we multiply (j2.15p by 
ii{e),e) foi" some function /3 G C^(M) with /3' bounded, and obtain 

^il^glA + B . V(,,^)/?(f ^(,),,) = :4e/3'(f ^(,),,). 
By letting e go to zero, we obtain the equation 

^ + ^.V(.,^)/3(f) = 0, 

in L°°([0, T]; n ^^^) for such functions /3. Now, letting (3 approximate the absolute 
value function, we end up with 

^ + i3.V(.,^)|f| = 0. 

This implies that we have a nonnegative solution |f | to (j2.ip . which vanishes at initial 
time and belongs to the functional space stated in Lemma |2.2[ Applying Lemma 12.21 we 
get |f I =0, that is, /i = f2- There remains now to prove the existence part. 

Existence in the functional space -L°° ([0, T]; L|^^ n L^^^(M^)) is given in a straight 
forward way by an application of Proposition 2.1 of [8j. For the sake of consistency, let 
us only mention here that it is a simple matter of regularization of the vector field B 
appearing in (j2.ip . That is, one introduces the solution fa to 

^ + ^„-V/a = 0, in (0,oo)xM^ 

where Ba = * B £ ^^([0, T]; C°°(IR^)) converges to B, then shows the desired estimates 
on /„, and finally passes to the limit. 

Next, the non-standard part we have to prove here is the fact that such a solution 
necessarily satisfies G L°°{[0,T], L^{R^). This is actually a consequence of the specific 
form of the transport equation and of the regularization process we have already done. 
Indeed, first, by the method of characteristics, we know if /o > a.e in M^, then f{t) > 
a.e in for all t > 0. Then, formally we multiply (j2.ip by to obtain 

Then we integrate the above identity over ^ on M'^ to deduce 
d_ 
di 

since 



/ \^ffdxd^= [ diY^m\E + ^xB))fdxd^, (2.16) 



[ lei'e • v,/dxde = - / diY,m^ofdxd^ = o. 

For the term on the right-hand side of (j2.16p . we have 

/ div^i\^\\E + ^xB))fdxd^ = 2 [ {^■Ef)dxdt 

since 

div^dePe xB) = 2C-{^xB) + |epdiv5(e xB) = 0. 
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Also, notice that, for a.e x G M^, 



m<R} Jm>R} 



<W3i?'||/||L-(M6)+i?-l 



m>R} 



(2.17) 



where ws is the vohime of the unit ball in M^, and in the last inequality R is taken to be 
Hence, we have the following estimate, by the Holder inequality. 



dw^m\E + ^xB))fdxd^ 



= 2 

< 2 

< C 

< C 



(C • Ef) dxdi 
\mE\dxd^ 

n^fdi] ' \E\dx 

4 

i\^fdxdA ||^||loo([o,t],l5 



Substituting this back to ()2.16p . we obtain 

^ f \C\'fdxdC<c( [ ICl'fdxdC] ||i?|[L°<>([0,T],L5(R3)). 



This implies 



[ \^\''f{t)dxdC< [ \i\^hdxdi + CT\ 



for all t G [0,r]. 

Finally, we show that the solution / necessarily belongs to L°°([0, T]; L!^(R^, L|(M^))) 
if /o G LI^(M^, L^(R^)). This is actually also a consequence of the specific form of the 
transport equation and of the regularization process as mentioned earlier. Indeed, we 
mollify E + ^ ■ B hy Ka obtain. 



dfa 

dt 



+ i ■ V,/, + {E + ixB)^- V^fa = 0. 



(2.18) 



Integrating (|2.18p over ^ in M^, one has, thanks to the fact that div^(£' + x B)a = 0, 

d 



dt 



fad^ + i-V^ / fo.di = Q 



I 



2.19) 
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That is equivalent to saying that fadi satisfies a conservation form, which yields the 
conservation over time of 



fadi 

Hence, G L°°{Q,T]Lf{W^,L^{^f))). By letting a ^ 0, one obtain 

/ G L°°(0,^;L^(M^Ll(M|))). 
The proof of Theorem 12.11 of complete. □ 



Remark 2.1. The assumption E G -L°°(0, T; L^(R'^)) is only needed to show the uniform 
estimate f G L'^{0,T;L^{ 



We now turn to the extension of the previous result to less regular initial data through 
the notion of renormalized solutions in the spirit of [8]. As in [8j, we consider the set 
of all measurable functions / on with value in M such that meas{\f \ > A} < oo, 
for all A > 0. For any /3 G Co,o(IK)) bounded and vanishing near zero, we thus have 
/3(/) £ n L°°{R^) for any / G L^. As in [8], we shall say that a sequence is bounded 
(respectively, converges) in whenever /3{fn) is bounded (respectively, converges) in 
for any such /3. But now we need some additional assumptions on our initial data, and 
that is why we consider the subset L^^ of consisting of functions / satisfying 

/ l^l"^ dxd^ < cs < oo, V<5>0. 

•''{|/(^,0I>'5} 

This subset is equipped with the topology induced by that of L^. For any / G L^^, we 
have |Cp/3(/) G L^^. ^^(M^). Indeed, for 6 small enough such that /? vanishes on [0,6], we 
have 



/ 



\c\''m{x,o)\dxdc= / \c\^mfix,mdxdc 

■J {\fi^,l;)\>S} 

+ / \c\''mix,o)\dxdc 

J{\f{^,0\<s} 
< ||/3||l-Q + < oo. 



It follows that if we choose /o in L^^, then /3(/o) is a convenient initial condition for the 
transport equation considered in Theorem 12.11 We therefore say that / is a renormalized 
solution of ()2.ip complemented by an initial condition /o G L^^ whenever /?(/) is a solution 
of ()2.ip in the sense of Theorem 12.11 with the initial condition /3(/o)- 



3. Stability of Vlasov-Maxwell-Boltzmann Equations: Main Results 

Let us begin by recalling that the general Vlasov-Maxwell-Boltzmann equations (jl.ip 
has the collison operator Q{f, f) which can be written as 

W,/) = Q+(/,/)-Q- (/,/), 

where 



Q^ifJ) 



[ dc [ dub{^-^,,u)f'f: 
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and 



Q-{fJ) 

with 



p3 



(3.1) 



L{f) = A*i:f, A{z)= / h{z,uj)du, 

The colhsion kernel b in the collision operator Q is a given function on x S'^ . We 
shall always assume the so-called angular cut-off kernel throughout the rest of the paper, 
that is, h satisfies 

b G L^{Br X S^) for all R G (0, 00), 5 > 
where Br = {z G : l^j < R], and 

{h[z,w) depends only on \z\ and |(2;,a;)[, 
(1 + /^^^^ A{^)di^ 0, as \z\ -> 00, for ah R G (0, 00). 

A classical example of such angular cut-off collision kernels is given by the so-called hard- 
spheres model where we have 

h{z,uj) = \{z,uj)\. 

The VMB system (jl.ip is complemented with the initial conditions 

f/|i=o = /o, on with /o>0, 

\E\t=o = Eo, B\t=o = Bo on Rl, 

with the usual compatibility condition 

divi?o = 0, and divE'o = po = / fod^, on M^. 

We state below our main stability results concerning the Cauchy problem of the Vlasov- 
Maxwell-Boltzmann system p.ip and p.ip . We assume that /o satisfies 

[ fo{l + iy+\^\^ + \logfo\)dxd^+ [ {\Eo\^ + \Bo\^)dx<^, (3.2) 

where v = v^x) is some function in satisfying 

1/ > 0, (1 + u)^ is Lipschitz on M^, e^" G L^{R^). 
Using the classical identity (see Lemma 2.1 in [3]), 

/ QU\mm = ] [ dxdiJ dojb{f'fi-ff,){c+c-c'-c), (3.3) 

we deduce the following local conservation laws of mass, momentum and kinetic energy: 

^ + dw,j = 0, (3.4) 

f^d^ + ExB) 

/ /Il2ll2 \\ V / 

+ dw^( [ g + ( 'l^'^^ Id-E®E-B(S)B)] =0, 
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^ iL ■^'^'''^^) ^ iL ^i^''^'^^) L ^-^"^^ " °' ^^'^^ 

for {x,t) G M'^ X (0,00). In fact, while p.4p and p.6p are easy to verify, we need to pay 
more attention to (j3.5p . To verify (j3.5p . we first multiply (jl.lap by ^ and integrate with 
respect to ^ to obtain 

I: / f^d^ + div, / e ® e/f^e = - (pi? + i X . (3.7) 

Note that 

EdivE + (V X £;) X = div(£; (g) E) - ^V|£;p. 
Thus it yields the following, combined with (jl.lbp and (jl.lcp . 

— {E xB) + div^ f ' ' ^ ' ' Id -E(gE-B®B\= -{pE + j x B). (3.8) 

Then, adding (I3.7p and (13.81) together gives (13. 5p . Integrating (I3.4p - (l3.6p in x over R^, we 
deduce the following global conservation of mass, momentum and total energy 

— / fdxdi = 0, for t > 0, (3.9) 

at 7r6 

T" ( / fidxd^, + f E X Bdx) = 0, (3.10) 

4- [ f\i?dxdi -2 I E- I ifd^dx = 0, for t > 0. (3.11) 

On the other hand, multiplying (jl.lbp by E, multiplying (jl.lcp by B, integrating them in 
X over M^, and then summing them together, we obtain 

d 



dt 



[ {\E\^ + \B\'^)dx = -2 [ E 



jdx. 



Substituting the above identity back to (j3.1ip . one obtains 

T" ( / f\i?dxd^+ f {\E\^ + \B\^)dx\ for t > 0. (3.12) 

Therefore, if we assume that the initial condition /o as (j3.2p . we deduce from (j3.9p . 
(j3li)]) and (IXT2D that 

sup / f{l + v + \i\^)dxdi+ [ QEl"^ + \B\'^)dx <C{T) (3.13) 

for some nonnegative constant C(T) that depends only on T and on the initial data. 
Indeed, we observe that we have, multiplying (jl.lap by ^{x) and then integrating over ^, 

^( [ fHx)dA + div, ( I fu{x)idi] = I fi- V,v{x)di 

< \ ! m?di+\p{t,x)\Vv\^ 
2 Jm.3 2 

< \ I m^di+c [ fd^+c [ fudi, 

2 Jr3 Jr3 Jk3 
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since {1 + 1^)2 is Lipschitz. In particular, we deduce 

^/ fu[x)dxdi<C + \ f fm^ + i^{^))dxdi. 

Then ()3.13p follows from the above inequality and Gronwall's inequality. 

The final formal bound we wish to obtain is deduced from the entropy identity. Multi- 
plying (jl.lap by log/, using (j3.3p . we obtain, at least formally, 

^/ flogfdxd^ + ] [ dx [ d^diJ Bif'f:-ff,)log^ = 0. (3.14) 

Since the second term is clearly nonnegative, we deduce in particular that 

sup/ f log fdxdC< fologfodxd^. (3.15) 
t>o Jr6 Jr6 

This inequality together with a lemma in [52] implies 

sup / /|iog/|<c(r). 

t€[0,T] JR6 

Also, if we go back to (j3.14p . we deduce that, 

f^dtf dx f didiJ 6(/7:-//*)log^<C(T). 
In conclusion, we obtain the following bounds: 



sup ( / f(l + \i\^ + v{x) + \\ogf\)dxdi+ [ (I^p + [B|2)dx) <C{T); 
t&[0,T] V^R6 JR3 / 

r [ dx f didiJ' dujbif'fi-ff,)\og^<C{T). 
JR» JR6 Js^ J J* 



(3.16) 



Now we give the definition of Renormalized Solutions to VMB. 

Definition 3.1. A triple (/(t, x, S^),E{t, x),B(t, x)) with / > is said to be a renormalized 
solution to VMB ([TT]) if for ah T e (0, 00), we have 

• f{t,x,C) G C{[0,T];L^R^)), E,B e C{[0,T]- L'^{R^)) and (l3J6]l holds; 

• for any /3 G C^([0,oo)) satisfying that /3(0) = and l3'{t){l + t) is bounded in 

[0,00), 

§-/{f) + e • V./3(/) + (E + ^xB). = P'UlQif, f) (3.17) 
holds in V (sense of distributions); and 

• (fTTb]) and (fTTc]) hold in V. 

One of the main objectives in the rest of this paper is devoted to the stability of 
renormalized solutions to VMB. More precisely, we consider a sequence of initial data 
{{fl^,El},BJ})}^^^ satisfying 1^ with f^^ > 0, a.e. in and converging to {fo,Eo,Bo). 
Then, corresponding to those initial conditions, we suppose that there is a sequence of 
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renormalized solutions {(/", -E", S")}^^ to VMB satisfying (j3.16|) . Without loss of gen- 
erality, we may assume that converges weakly to {f,E,B). We will prove 
that {f,E,B) is still a renormalized solution to VMB with the initial data {fo, Eq, Bq). 

Theorem 3.1 (Weak Stability). Suppose that {(/", i^'", -B"')}5^i is a sequence of renor- 
malized solutions to VMB (jl.ip satisfying (j3.16p . with initial data {(/q j Eq, Bq)}^^^ satis- 
fying (I32D, fo > 0, a.e. in and converging weakly to (/o, £^0, Bq) in L^{M.^) x (L'^iM.^))^; 
and{f,E,B) is a weak-* limit of {{f' , E"" , B'')} in L'^{0,T; L\R^))x (^L°°{0,T; (L2(r3))^^. 
Then the sequence {/«} satisfies: 

(1) For all ip G C(M3) such that — > as |^| — > oo, j^j, f^cpdS, converges to 

^3 /V'de in LP{[0,T],Ll^{R^)) for alll<p<oc. 

(2) L(/") con?;er5es to L(/) m LP([0, T]; L1(M3. x K)) for alll <p < oo, T £ (0,00), 

compact set in M|. 

(3) For all (j) G L°°(M3) wit/i compact support, J^^ Q^{f"', f"')4>d$, converges locally 
in measure to J^s Q^{f, f)4'dC- ^nd Q^{f"', f"'){l + f"')~^ are relatively weakly 
compact in L^(M^ x K x (0,T)) for all T G (0, 00), compact set K in M|. 

(4) Q~^{f^ , f^) converges locally in measure to Q^{f,f)- 
Moreover, if 

II^"IIl°°(o,T;L5(r3)) is uniformly bounded, 

and 

ll-^"llL°°(o,r;L''(R3)) is uniformly bounded for some s > 5, 
then the weak limit (/, E, B) is a renormalized solution of (jl.ip with the initial data 
{fo, Eq, Bq). 

Remark 3.1. Due to the convexity of xlnx and the monotonicity of (x — y) In | for all 
x,y > 0, we can show, as in [llj . 

I f{t)lnf{t)dxd^<limmf [ /"(t) In /"(t)dxdC, 

and 

f f [ didiJ dubif'fi- ff,) In ^ 

<liminf/ / / didiJ dub{f^'fK-rf:)ln^-f-^, 

for all t > 0. This entropy estimate is crucial for the long time behavior of renormalized 
solutions. 

A consequence of the weak stability is the propagation of smoothness of renormalized 
solutions. 

Theorem 3.2 (Propagation of Smoothness). If, in addition to the assumptions in Theo- 
rem \3.1i f^ converges in L^iR^) to fo, then /" converges to f in C([0, T]; L^(M^)) for all 
T G [0,00), and {f,E,B) is a renormalized solution of (11. ip if {f^ , E^ , B'^) is a sequence 
of renormalized solutions. 
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Remark 3.2. The assumption that E{x,t) is uniformly bounded in L°°(0, T; L^(M'^)) is 
crucial for Theorems l3.1l and l3.2l because of the nonlinear term associated with the Lorentz 
force. Notice that usually from Maxwell's equations, we can only obtain the a priori 
estimates on E and B in L'^{0,T; L^{R^)). 



4. Proof of Theorem I3.lt Weak Stability 

This section is devoted to the proof of Theorem 13. II We divide the proof into two steps. 
In the first step we show why the first four statements of the theorem hold. Then we 
concentrate in the second step on the proof of the fact that the weak limit is indeed a 
renormalized solution of Vlasov-Maxwell-Boltzmann equations. We remark that the first 
step is essentially an adaptation of the results and methods of [lOl ESJ [23] , while the second 
one requires a new result of renormalized solutions for the Vlasov-Maxwell equations. 

4.1. Step One. In this subsection, we are aiming at proving the first statement of The- 
orem 13.11 following the spirit of [10] . Then the second and the third statements can be 
shown exactly as in [10]. Finally, once the first three statements hold, the fourth statement 
will immediately follows from the argument in [22]. Therefore, for the sake of conciseness, 
we only give the detailed proof of the first statement of Theorem 13.11 

In order to prove the first statement, we first recall that for all compact sets C M| 
and T G (0, oo), we have 

/ (1 + n-'Q-ir,ndxd^ < [ Lindxd^ 

= [ dx [ r(x,e*,t) / (4.1) 

<C [ dx f r{x,C„t){l + |e*P)de* < oo, 

due to the assumption on the collision kernel b, hence, 

(1 + ry^Q'{r,n is bounded in L°°{0,T;L\R^ X K)). 
Also, we observe that we have (see [TO]). 



Q+ir, n < 2Q-(r , n + T^i d^* [ u^ir' - /: - ri:) in 

in 2 Jk3 Js2 



which, combining with (j3.16p and (j4.ip . implies 

{1 + fT^Q^{r,r) isboundedin L^(0, T; ^^(M^ x i^)) (4.2) 

for all compact sets K in ]R| and T € (0, oo). 

Next, we observe that since /" is a renormalized solution of VBM (|1.1|) . we have, for 

^ + e • V.) = f^sifiQir, n - div^((s" + e x B^)Ps{n) (4.3) 

in v. In order to apply the velocity averaging results in [Oj |T3], we remark that ()4.ip 
and (jMl) imply that Psir)Q{r,D is bounded in L^{0,T; L^{M.l x K)) for all compact 
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subsets K of M|. And also we observe that l3s{D is bounded in L°°{{0,T) x R^)), and 
hence, div^((£;" +^ x B'^)f3siD) is bounded in L^{iO,T) x M^; i7-i(M3)). Denoting 

rsin = miQirjn, 

and decomposing f3s{f^) into 

= Mnx{(t,.mrsin\<M} G L^mr) X M.')), 

g"-, and /i" by 

+ ^ ■ ^0 ='7'5(r)X{{t,.,0:|r,{/")|<A/} ^^^^^ 

- div^ ((i?" + C X i?")/35(r)X{(M,?):|r.(/")|<M}) , 

+ e • V,) <7" = -div5((i?" + e X i3")/35(r)x{(M,5):ir.(/")i>Af}), (4.5) 

(I + e • V.) = Ts{nx{it,.,mrsir)\>M}, (4.6) 

for M > 1, where 

/i"|t=o = 5"k=o = 0, n"|t=o = /35(/o), 
and X is the characteristic function of sets. Because {Ts{f^)}'^=i is weakly compact in 
L^((0,r) X R^) due to the facts that /3'(t) = jj:^ < and j:^Qir, P) is weakly 
compact in L^((0,T) x M^), and because, from ()4.6p . 

+ = / 75(/")x{(t,x,0:|ri(/")|>M}(T,2; + ^r,Odr, 

JO 

it follows that, uniformly with respect to n, 

I ff \hJ'{t,x,i)\d^dxdt ^{), as M ^ oo. (4.7) 

Similarly, from the compactness of Ts{f^''), we deduce that 

5- := (ii;- + e X i?")/35(r )x{(t,.,0:|r.(r)l>M} ^ o (4.8) 
in L}^^{{0,T) X R^) as M ^ cxd. From (l45|), we have 

X + te, = f -div5((i?" + e X i3")/35(r )x{(i.x.O:|r,(/")l>M})(r, x + ^r, O^r. 

JO 

Thus, for any -0 S Pg(M'^), we deduce from the above identity that 

= f I {{E^ + e X B^)Ps{nX{it,xmrsU-)\>M}){r,x + er,0 • V^^d^dr. 
Jo JR^ 

Therefore, from the weak compactness of 5", the above identity with (j4.8p implies 



in LU(0,T)x 



/ c/"(t,x,OV'c^C ^ 0, as M^oo, (4.9) 

JR3 
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On the other hand, since {u"'}'^^^ and {Ts{f^)x{{t,x,£,):\Ts(f")\<M}}'^-i are bounded 
sequences in L'^{{0,T) x R^), and div^{{E'' + ^ x B"')l3s{r)) is bounded in L'^{{0,T) x 
R3;i/-^(R3)), by the velocity averaging lemma (Theorem 3 in [5]), we deduce that 



/ 



u"-ip{Od^ is bounded in H4{{0,T) 



p3n 



for all V G V{R^). Thus, {/jj3 w"^(Oc^C}r=i compact in L'^{{0,T) x R^) and is locally 
compact in L^((0,T) x R^), which, combining with (14. 7p and (I4.9p . implies that 

/ P&iDi^di is relatively compact in LP(x (0, T), ^/^^(M^)) (4.10) 

for ah 1 < p < oo, G D(R3). 

The first statement of the theorem for -0 € P(R^) then follows from (|4.10|) and (|3.16|) . 
since it suffices to observe that we have for all i? > 1, 

< r - Psin < sRr + rx{r>R} < snr + r (4.11) 

and then take the limit as — > oo and 6^0. Next, for a general € C(R^) such that 
+ ICP)^^ as 1^1 —7- oo, we introduce 



VM = 1] 



for M > 1, where t] G D(R'^), 0<??<1, ?? = lon Bi. Then the first statement holds for 
iprjM, and the first statement will be valid for such a tjj provided that 

sup/" dt [ - rjM)dxd^ ^ as M ^ oo, (4.12) 

n Jo JkxR^ 

for compact subsets € Ri^. Indeed, (j4.12p follows from p.l6p since 

nm-mi)dxdc 



JKxR-^ 

T 



< c sup [ dt [ r (1 + iep)x{ici>M} 



_\m\ 

l5l>A-f 



for some C > independent of n. 

4.2. Step Two. We now aim at proving that (/, E, B) is a renormalized solution of VBM. 
First of all, we claim that it is enough to show that 

Lemma 4.1. /// G L°°(0, T; L1(R6)), the equation (pDTD holds if and only if 

^ ln(l + /) + div,(C ln(l + f)) + {E + ^xB)-V^ ln(l + /) = J^^QU, f), (4-13) 
in v. 



WEAK STABILITY FOR THE VLASOV-MAXWELL-BOLTZMANN EQUATIONS 21 

Proof. On one hand, if / is a renormalized solution to VMB, then (j4.13p automatically 
holds since /?(/) = ln(l + /) G C\[0, oo)) with /3(0) = and (3'{f){l + /) = !. 

On the other hand, if ()4.13p holds, we claim that / is a renormalized solution to VMB. 
Indeed, denoting 

a{s) = /3(e^ - 1) 

for aU (3{t) € C^{[0, oo)) with /3(0) = and /3'(t)(l + t)<C. Then, we have 

dMf) = a'{f)dtf; 

Vxa(/) = a'(/)V,/; 

V«a(/) = a'(/)Vg/; 
Multiplying dH]) by cr'(ln(l + /)), we obtain, 

^a(ln(l + /)) + div,(ea(ln(l + /))) + {E + ^ x B) ■ + /)) 

= a'(ln(l + /))^g(/,/), 

in the sense of distributions. Note that, by the definition a, we have 

a(ln(l + /))=/3(/), and cT'(ln(l + /)) = /3'(/)(l + /). 
Hence, substituting the above two identities in (|4.14p . we get 

^/3(/) + div,(e/3(/)) + {E + ixB)- V^/3(/) = /3'(/)Q(/, /), 
in the sense of distributions. The proof of this lemma is complete. □ 

The rest of this subsection is devoted to the proof of ()4.13p . Recall that we deduce, 
from a priori estimate (I3.16P and weak passages to the limit. 



(4.14) 



sup ([ f{l + \i\^ + v{x) + \\ogf\)dxdi+ [ {\E\^ + \B\ 
te[o,T] \jk6 Jm? 



')dx 



+ r f dx f / d^6(/7:-//*)log44 <oo. 



(4.15) 



for all T e (0, oo). Now the strategy to prove ()4.13p is the following: we first consider 

/35(r) = r(i + <5n-' 

for 5 € (0, 1] and weakly pass to the limit as n goes to oo in the equation satisfied by 
Psif^)'^ then for the equation satisfied by the limit of Psif^) as n — ?> oo, we use /3 to 
renormalize it and let 5 go to to recover (|4.13p . To begin with, without loss of generality, 
in view of (j3.16p . we can assume 

weakly* in L°°(0, T; L^(M^)); 
B'^^B weakly* in L°°(0, T; ^^(M^)). 
E^^E weakly* in L"^ {<d,T-L^{W^) r\ L^{W^)). 

Furthermore, without loss of generality, extracting subsequence if necessary, we may as- 
sume that for all 5 > 

P&in^Ps weakly in LP(]R6 x (0, T)); (4.16) 
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/i^ = (1 + sry^ hs weakly-* in L^{M.^ x (0, oo)); (4.17) 
g- = r{l + 5n-^^gs weakly in L^iR^ x {0,T)), (4.18) 

for all T £ (0, oo), 1 < p < oo. Furthermore, because of the third statement and the 
equi-integrability, we may assume that 

{l + 5n~^Q^{r,n^Qf weaklyln ^^(R^ x x (0, T)), (4.19) 

for all compact sets C M| and T € (0, oo). 

Notice that, since /" is a renormalized solution of VMB, (|4.3p holds with /3(/") replaced 
by I3s{f^) for all 5 > and we want to pass to the limit in these equations as n goes to oo. 
To this end, we deduce from the first statement of Theorem 13.11 that and j" converge 
in LP(0,T; L^(M^)) to p and j, respectively for all 1 < p < oo and T E (0,oo). We then 
pass to the limit in (|4.3p and we obtain 

^/35 + div^(e/35) + (^ + exS)-V/35 = Q+-Q7, xGM^ ^ G M^ t > 0, (4.20a) 
dE 

_ - V X S = -j, divB = 0, on M^. x (0, oo), (4.20b) 
dB 

— + V X E = 0, dwE = p, on M^. x (0, oo), (4.20c) 



P 



I fdt 3=1 fidi, on m3x(0,oo), (4.20d) 

in T>' . Here, for the convergence of the nonlinear term (i?" + ^ x i?") • V/35(/"'), we need 
to show, for all (f) E ^^((0, oo) x 

rt /> 

(^(^" + ^ X B") • V^/3sir)dCdxds 



oJu^ (4.21) 
/ Vgt/) • (^" + e X B'')^s(.r)d^dxds, 

since 

div5(^" + e X 5") = 0. 

If we take (j) = 0(t,x)<I>(^) (which is enough by dense property) for G P((0,oo) x M.^) 
and <I> G P(M'^), we can rewrite the term on the right-hand side of (|4.21|) as 

- f f Ht,x){E^+^x B'^) ■ ( f mMfldA dxds, 
Jo Jr3 Vjm3 / 

by letting tp = V^<I>. In fact, on one hand, by ()4.10p or the velocity averaging lemma 
in k^HOMndC and J^^mOMndC strongly converge to f^sHOMC and 

/iR3 Ci^iOf^sdC in LP(0, T; Lj^^(R^)) respectively. On the other hand, since G Co{R^) and 
f3s{t) < t, we have, using (|2.17p . 



/ Ci^(3s{r)d^\ is uniformly bounded in L~(0, T; Li (M^)), 

J n=l 

and ^ 

{/ V'/?5(/")dQ is uniformly bounded in L°°(0, T; lI (M^)). 

UM3 J „=i 
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The latter is true, because, for all i? > 1, 



where |-B(0, 1)| denotes the Lebesgue measure of the unit ball -6(0, 1) in M^, and by taking 

's\\\er 



R 



|2 fnll \ ^ 



(|4.22p becomes 



5(0,1)1 ; ' 

[\B{oA)\h\\enlu^^ 



\ 



5l 



Therefore, 



%{r)di \ is uniformly bounded in L°°(0, T; (R^)), 

J n=l 

since {|CP/"}^=i is uniformly bounded in L°°(0, T; L1(M6)). Thus 

/ i^^sindi^ 1 i-^Mi in LP(0,r;Lf„,(M3)) for all 1 < s < ^, 

and 

/ ^PsiDdi^ f i^Mi in LP(0,T;LL,(M3)) for all l<r<\, 
for all 1 < p < oo. 

The weak convergence of in L^((0,T) x M^), combined with the strong convergence 
of iRsi'iOl^siDdt implies 

/ / V /:<!)■ E"^ PsiDdidxds ^ [ [ V /:<!)■ Epsd^dxds. 

Jo JR6 Jo Jr6 

The similar argument goes to the second part of the nonlinear term 

^{t,x)B^x( [ iPiOmnd^) dxds, 



due to the weak convergence of in L^((0,T) x M'^) for q > 5. That is, 

[ [ V5</> • C X B'^psiDdidxds [ [ V^(l)-^x BPsd^. 
Jo Jr^ Jo Jr^ 

Next, since /35(/") € ^^(RS) nL°°(M6), we know that e (M.^) D (R^) . Also, since 

lePr e L-(0,r;Li(R6)), we know that PsiDlCl^ G L-(0, T; L^M')) and {Mr)}n=i 
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is weakly compact in L°°{0,T; L'^{R^)). Hence \C\'^f3s e L°°{0,T; L'^{R^)). Thus, for any 
(7 > 0, we have 

/ jel^dxd^ < - / l^l'^Mxd^ < - [ iCfPsdxd^ < oo. 

Therefore, Theorem 12.11 impUes that 13s is a renormahzed solution of (I4.20p . 
As (5 — )• 0, we claim that 

Lemma 4.2. 

(3s ^f, in C{[0,T];L\R^)), 

as 6 ^ 0. 

Proof. We start with proving the continuity of f3s with respect to t > with values in 
LP{M.^) for all 1 < p < oo. To this end, we remark that if we regularize by convolution I3s 
into /3| as in Lemma |2.H we obtain 

^/3| + e • V,/3| + {E + CxB). = Qj- QJ + (4.23) 

where ^ in L^{0,T; L}^^{R^) n L°°(M'5)) as e goes to for ah T G (0,oo). Hence, it 
is easy to see from (03]) that, /3| G C([0, 00); LP(M6)) for 1 < p < 00. Note that (3s is a 
renormalized solution to the VM (|4.20ap . Subtracting (j4.23p from (j4.20p . multiplying the 
result by \(3s — (3^\p~'^{(3s — and then integrating over M^, we obtain 

-T \^s - ^ll^'dxdC ^ in L^{0,T), as e^O (4.24) 
dt J^6 



for all 1 < p < 00, r G (0, 00). It follows that ^s e ^([0, T]; L^i 

Next, we show that / G C([0, 00); L^(M^)). Indeed, because of (j3.16p . we have for all 
T G (0,00), as in 

sup sup|[r-/35(r)||ii(M6) ^0 as 5^0. (4.25) 

te[0,T] n>l 

Hence, by the lower semi-continuity of the weak convergence, we obtain 
sup \\f - (3s\\L^(Re) < sup liminf ||/"-/35(/")||ii(iR6) 

< sup sup||r-/35(r)||Li(R6)^0 as 5^0, 

tG[0,T] n>l 

and this implies (3s converges in C([0, T]; L^(M^)) to /. □ 

Now we can state the equation (|4.2Uap more precisely. To this end, we observe that 
TT3t' {i+St)"^ ^'■^ convex on [0,oo), therefore we have 

< (^sif), hs>{l + Sf)-"^ a.e on R^ x (0, 00). (4.26) 

In addition ^i_^^^yi = (3s{'t)i^ — S(35{t)), because the function x{l — 6x) is a concave function, 
hence 

9S < Ps{l - S(3s) a.e on R^ x (0, 00). (4.27) 
Furthermore, because of the second statement of Theorem 13. H we deduce that 

Qj=gsL{f) a.e on M^x(0,oo). (4.28) 
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And, using the fourth statement of Theorem 13. H we could also deduce that 

Qj = hsQ+{fJ) a.eon R'^x(0,oo). (4.29) 
We finally use the fact that /J^ is a renormalized solution of (|4.2U|) to write 

= {l + l3s)-'Qj-{l + f3s)-'Qj\ 

And we wish to recover (j4.13p by letting 5 go to 0. Recall that we already showed in 
Lemma [O] that converges to / in C{[0,T]; (M.^)) for ah T G (0,cx)). Therefore, in 
order to complete the proof of Theorem 13. 1^ it only remains to show 

Lemma 4.3. 

Qf{l + f^s)^'^ are weakly relatively compact in L^{M.l x K x (0, T)) (4.31) 
for all compact sets K CM^^ and T G (0, oo), and 

{l + (3sr^Qj ^{l + f)-^Q-{f,f), a.e (4.32a) 

(1 + f3s)-'Qj ^ (1 + fr'Q^if, /), a.e (4.32b) 

as 6 goes to 0. 

Proof. We will follow the lines of the argument in [23] and begin with Qj . Without loss 

of generality, we may assume that /3s converges a.e. to f as 6 goes to 0. Then, (|4.32ap 
follows since 

(1 + i^sr'Qj = (1 + Psr'gsLif) ^ (1 + fr'fLif) 

(5 — 7> provided we show that gs converges a.e. to /. 
This is easy since we have for all > 1, 

< r - r (1 + sf^-' < sR6r + rx{r>fi}, 

hence gs converges to / in C([0, T]; L^(IR^)) for all T € (0, oo) by the uniform integrability 
of /" and the lower semi-continuity of the weak convergence. We now prove ()4.3ip for Qj 
by first observing that (|4.28p yields 

< (1 + Psr'Qj = (1 + M'^gsHf) 

<i^-mYf^m<L{f), a.e. 

And we conclude the proof of (|4.3ip for Qj by the equi- integrability, since L{f) G 
L°°(0, T; L^{Rl x K)) for all compact sets iiT C and T G (0, oo). 

Next, we turn to the proof of ()i:3T1) for and ()4.32bll . We begin with (031]). We 
recall the following classical inequality for all M > 1, 

g+(r,r) < MQ-(r,r) + (4.33) 

where 



is positive and bounded in L^{M.^ x (0,T)) for all T G (0,oo). Without loss of generality, 
we may assume that e" converges weakly in the sense of measures to some bounded 
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nonnegative measure e on x [0, oo) and we denote by cq its regular part with respect 
to the usual Lebesgue measure, that is, cq = j^yh' iu^^) £ x (0)^)- Dividing (I4.33P 
by (1 + (5/")^ and letting n go to oo, we obtain 

hence 

Qj <MQJ + j^eo a.e. on W x (0, oo). 

Then (14.310 for Q'^ follows since we already show it for Qj and the integrability of bq. 
We finally prove ()4.32p for Q'^ . We first remark that we have for all i? > 0, 

Q+(r , n > (1 + sr^Q^ir, n 

o _u (4-34) 

>(.i + 5Rr'Q+{r,nx{f-<R}. 

In particular, if we multiply (j4.34p hy ip & Cq°(M2) with ^ > 0, we find by letting n go to 
oo and using the third statement of Theorem 13.11 

/ Q^ifJWC> [ Qj^di a.e. on ^^^(Coo). 

Indeed, the integrated left-hand side converges locally in measure while the right-hand 
side converges weakly in and this is enough to pass to the limit in the inequality a.e. 
on X (0, oo). Therefore, we have for all 5 € (0, 1], 

Q+(/,/)>Q|. (4.35) 

Next, we use the other part of the inequality (j4.34p and we write for r G (0, 1], using 
(1133]), 

{i+5R)~\i+TL{n)-^Q^u\n 

< (1 + 5r )-2Q+(r , n + (i + rL(r ))- V{/">i?}Q+(r , n (4.36) 

1 „ M 

e -I 

InM T 

We then observe that Q^{f "', /")(1 + TL{f'^))~^ is relatively weakly compact in L-'^(M^ x 

InM* 

we may assume without loss of generality that it converges weakly in L^(M^ x (0, T)) for 
all T G (0,00). We claim that its weak limit is given by (1 + TL{f))^^Q^{f, f). Indeed, 
ii ip & L°°(IR|) with compact support, we have 

(i+rL(r))-iQ+(r,r)V'de= / g+(r,r)vx, 



< (1 + srr'Q^r, n + j^e- + ^/-xir >«>. 

^rve that /")(1 + rL(/"))-i is relatively weak 

(0, T)) for ah T G (0, 00) since it is bounded above by j-^e^+MrP for all M > 1. Hence, 



where i/jj^ is uniformly bounded in L°°(IR^), and has a uniform compact support and 
^tP^ = (l + rL(/))~V in LP{{0,T) x M^) for ah 1 < p < 00. This is enough to enable 
us to deduce 

locally in measure on x [0, 00), which yields the claim. 
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We then pass to the hmit in (j4.36p and deduce as above 

(1 + 6Rr\l + TL{f)r'Q+{f, f)<Qt + + ^fn, a.e., (4.37) 



/ fRdxd^= lim / PX{f">R}dxd^ < 



where fn is the weak limit of f "'X{f">R} iii i^(IR^). Since we have 

C 
h^' 

we deduce from (j4.37p . by letting first J go to 0, then M go to c«, then go to oo, and 
finally r go to 0, that 

g+(/,/)<limQ+ a.e., 
which, combining with (|4.35p . imphes that 

Q+(/,/) = limQ+ a.e. 

o— >0 

The proof is complete. □ 

Putting together the conclusion of Step One, Lemmas 14. 1114. 3^ we finish the proof of 
Theorem 13.11 



5. Proof of Theorem 13. 2t Propagation of Smoothness 

In this section, we prove Theorem l3.2l First, without loss of generality, in view of (I3.16p . 
we can assume 

weakly* in L^{0,T; L^R^)); 
B'^^B weakly* in L~(0, T; l2(m3) n L^(M^)); 
E^^E weakly* in 1"=^ {<d,T-L^{W^) r\ L^{W^)). 

Applying Theorem 13. H we know that / € C([0, oo); L^(M^)) is a renormalized solution of 
VBM. In particular, we know that we have, setting 75(/) = \ ln(l + 5f), 

= i5{f)Q^{f,f)-Himf), 

in V'. It is easy to deduce that, 75(/) € C([0, oo); LP(M^)) for all 1 < p < oo since 
75(/) G C([0, oo); L1(M6)) n L°°(0, oo; L\r')), 

hence 

75{f)\t=o = isifo) a.e. on 
The strategy of the proof of Theorem 13.21 goes as follows. First of all, we introduce, 
without loss of generality, the weak limit of 7<5(/"') in L^(M^ x (0, T)) for all T € (0,oo) 
and 1 < p < oo, and we denote it by 75 (note the difference from the notation ■ysif"') 
throughout this section). The first step is to show that 75 is a supersolution of (15. 111 . 
In the second step, we deduce that 75 = 75(7) and that /" converges to / a.e. or in 
L"'^(M® X (0, T)) for all T € (0,oo). Finally in the third step, we show that /" converges 
to / in C{[0,T];L^{R^)), thus proving Theorem^ 
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Applying Theorem 13 . 1 1 and a similar argument in Section 4, we can show that 75 satisfies: 
75 G L°°(0, T; LP{R^)) for ah T G (0, 00) and 1 < p < 00, 

0<75<75(/) a.e. on m6x(0,oo), (5.2) 

and 

^ + div,(e75) + {E + ^xB)- Vas = Qt- QJ, (5.3) 

in D', where Q'^ , Qj are respectively the weak limits in L^(M^ x x (0, T)) for all compact 
sets C R| of (l + -5/")-iQ+(/", /"), (l+5/")~^Q"(/", /"). For the weak limit function 
75, we claim 

Lemma 5.1. 

75GC([0,oo);LP(M6)) 

for all 1 < p < 00. 

Proof. In fact, we claim that the weak limit 75 is a renormalized solution of (|5.3p and then 

75 gC([0,oo);LP(M6)) 
for all 1 < p < 00. For this purpose, we introduce 

7Kn = 7smfn) 

for e G (0, 1] and denote its weak limit by 7|. Then, the proof in Section 4 applies and 
shows that the weak limit 7I G C([0, 00); L^(M^)) is a renormalized solution of 



-7| + div,(C7|) + (i? + e X i?) • Vnl 



5_ 

= imifmif)Qt.e-^5i(^e{f)mf)Qj,e, 



where the notation g means the weak limit of the sequence {gn}^=i in ^joc- Next, we 
claim 

o<75(r)-7i(r)<r-/3e(n^o in lhr^) (5.5) 

uniformly in n > 1, t G [0,T]. Indeed, since the sequence {fn}^=i is equi-integrable, for 
any f] > 0, there exists two positive numbers D and R such that 



and 



sup / f^dtdxdS, < r], 

neN J([0,T]xBrxBr)'' 



sup / 



f'dtdxd^ < rj. 

{/">£»} 



Hence, in particular. 



sup / f^dtdxdS, < rj. 

n€N J{f">D}n[0,T]xBiiXBii 
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Therefore, we have 



sup/" ir - Pe{r))dtdxdi 

< sup / pdtdxd£, 

n&N J(\0,T]xBnxBT,Y 



n€N J{[0,T]xBiiXBiiY 

+ sup / rdtdxd^ (^-^^ 

neAf J {f">D}n[0,T]xBRxBR 

+ sup / ir - f3e{r))dtdxd^ 

neN J {fn<D}n[0,T]xBRxBR 

<2ri + D'^R^Te, 

since /" — /3e(/") < D^e if /" < D. Thus, letting first e go to and then r/ go to in 
(|5.6p . we deduce (|5.5p . 
Similarly, we have 

< 1 ^ in L1(M'') 
uniformly in n > 1, t € [0, T]; 

< (7K/5.(r)) -7Kr))/3;(r)Q-(r,r) 

Li(0,r;Li(M3 xK)) 
uniformly in n > 1, for all compact sets IT C M|; and 

< Hif^sin - isinWeirm^ir, n 

<Y^^Q-'{r,n^O in L\0,T;L\RlxK)) 
uniformly in n > 1, for all compact sets C M|. Here, we used 

Thus, letting e go to in (j5.4p . we deduce that 75 is a renormalized solution to ()5.3p . 



Hence, from (j5.3p . we deduce that 

^ G Li(0,r;Ty-"'i(M3)) 

for n > large enough. Also, we know that, since 75 (t) is a strictly concave function, 

< 7<5 < 75(/) < / e L-(0,r;Li(M6)). 

Hence, by the Aubin-Lions lemma in [21j . we know that 

-fseC{[0,T];W-''\R^)). 

But actually, we know 

75GL°°([0,r],L?'(M6)) 
for all 1 < p < cxo. Thus, by the interpolation, we know that 

75eC([0,T];L^'(M6)) 
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for all 1 < p < oo. □ 

5.1. Step One: 75 is a supersolution of (|5.1|) . Without loss of generality, we may 
assume that we have 

^s(f'') = YTJF ^ weakly* in L-(m6 x (0, 00)), 

and 

f'^^sif'") = YTJF ^ weakly* in L'^{R^ x (0,^)). 

Furthermore, since 7^(/), —t^'gif) are convex on [0, 00), we deduce the following inequal- 
ities: 

Cs>Y^=^s(.fl Os<-^ = H{f), (5.7a) 

75 < 7 ln(l + df) = -fsif) a.e. in x (0, 00). (5.7b) 


We claim 
Lemma 5.2. 

Q- = esL{f) a.e. on ]R^x(0,oo). (5.8) 

Proof. In fact, it is enough to verify that (|5.8p holds in [0, T] x Bji x Br, where Br is 
the ball with radius R and centered at the origin in M^. Due to the second statement of 
Theorem l3.lt we know that L(/") converges a.e. to L{f) in [0, T] x Br x Br. By Egorov's 
Theorem ([25]), for any e > 0, there exists a subset E C [0,r] x Br Br with \E\ < e 
such that L{f^) converges uniformly to L(/) on E^. Thus, for ah (j) G L°°(M^ x (0, T)), 

< Uh- sup f |L(/")| + \Ss\\L(f )\dxdidt 

n J E 



+ M 



Os L{f)dxdids 



+ \mL^\E'\sML{n-L{f)\- 

E<^ 

The first term can be made arbitrarily small uniformly in n, due to the equi-integrability 
of {L(/")}~^p The second term also goes to since L(/) G L^((0,r) x Br x Br). And 
the third term goes to as n goes to 00 since the uniform convergence of L{f^) to L{f) 
in E'^. Thus, (15. Sp is verified. □ 



Similarly, we have 
Lemma 5.3. 

Qj = C5Q^{f,f) a.e. on ^^^(Coo). (5.9) 

Proof. Indeed, let A be an arbitrary compact subset of x [0, 00). By the Egorov's 
theorem and the fourth statement of Theorem 13.11 for each e > there exists a measurable 
set E with the measure of E not greater than e (i.e., \E\ < e), up to a subsequence 
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Q~^{f^,f^) converges uniformly to Q'^if,/) on and Q^{f,f) is integrable on E'^. 
Then, for all (p G L°°(M^ x (0, oo)) supported in A, we have 

HisifiQ^irji - csQ-'u.mxdidt 

A 

<\mL^ I Ws{nQ^{r,n-csQ^{fj)\dxd^dt 



+ 



E 



<PW5in-Cs}Q-'{f,f)dxdCdt 
E<^nA 



+ Uh^ \E' n ^1 sup |Q+(r , n - Q^if, f)\, 

where the third term goes to as n goes to oo, for each e > by the uniform convergence 
of Q^U^i /") to Q~^{f-, /) on £"^. And so does the second term since 

'/'Xi?=Q+(/,/)G^'(K'x(0,oo)). 

Finally, since 'ys{f")Q~^{f^,f") is weakly relatively compact in L^(M^. x K x (0,T)) for 
all compact sets iiT C M|, the first term can be made arbitrarily small uniformly in n if 
we let e go to 0. 

Notice also that CsQ~^{f,f) S -^^^(IKx x K x (0, T)) by following the similar argument 
as before, we can show that C5iQ~^{fi f) A R) is tbe weak limit of i_^_]j„ /"") A R), 

where a Ab = min{a, b}. Thus, (j5.9p follows. □ 

Now, we use (ETD-fEni) in (O) to obtain 



^ + dw^iCls) + {E + ^xB).Vjs> l5{f)Q{f, f) (5.10) 
in V . We conclude this first step by proving that 75 satisfies the initial condition: 

75|t=o = 7<5(/o)- 

Indeed, in view of the equation satisfied by jsif^), we know that 

^^^Qp e Li(0,r,VF-"'i(M6)) 

for n > large enough, which, combined with the fact jsif^) e L°°(0, T; L^(]R6)) and the 
Aubin-Lions lemma, implies that 

75(r)^75 in Ci[0,T];W-''\R^)) 

for any s > 1. But, by the assumption, 75(/")|j=o = isifo) converges in L^(M^) and thus 
in W~^'^{M.^) to 75 (/o)- Thus, we conclude that 75 satisfies the initial condition. 

5.2. Step Two: 75 = 7,5 (/) and /" converges in to /. To this end, we consider 

75(/)-75 = r5eC([0,oo);LP(IR6)) 
and observe that ts satisfies, in view of (jS.ip and (jS.lOp . 

+ div,(er5) + {E + ^xB)-Vi:T5<0 (5.11) 
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in T>' with 

T(5 > a.e. on x (0,oo), Ts\t=Q = a.e. on 
Then, for ts we have 

Lemma 5.4. 

r5 = 0. 

Proof. Formally, we only need to integrate ()5.1ip over to get 



(5.12) 



d f 

— Tsdxd^<0 in ^'(0,00). 
at J^a 



(5.13) 



Then (f5l^ with (f02]) yield: T5 = on x (0, 00). 

Our main objective now is to justify ()5.13p . In order to do so, we introduce the function 
with 

Hz) 



'1, if \z\ < 1; 
0, if 1^1 > 2. 

Notice that /^^(t^) = also satisfies ([511]) and (I5l2]l . and /^^(r^) G C([0,oo);LP( 

for 1 < j5 < 00 and e > 0, since 

mx)-Pe{y)\ <\x-y\ 

for all x,y > 0. Then we multiply (j5.1ip by integrate the resulting 

inequality over x (0, t) for all t > to obtain 



/ I3s{ts)cP (-) cP ( ^) dxd^ 
< [ ds [ dxdC/3,iTs)-(^-V^(- 



(5.14) 



n 



- H- 

n \n 



Recall that 



sup (3e{TsM\'^dxdC : te[0,T], e > 0, 5 > o| 



< sup <j / /ICpfixd^ } < 00, 

tG[0,T] 



(5.15) 



since /Se^Ts) < ts < Jsif) ^ / for all T G (0, 00). Hence, for the terms on the right hand 
side of (I5.14p . we have 

ds / dxd^l3e{Ts) 



n \n 



n 



< 



< 2 



JR^ 



(t5 ) X{n< Ix 1 <2n} 2 1 1 1 1 ! I V0 1 1 Lex, dxdC 
fX{n<\x\<2n}'^\\4>\\L^ \\V dxd^ 
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as n — >■ oo and 

ft 



{E + ixB)-V(t)[- 

n 



[ ds [ dxd^l3e{Ts)-(l)(- 

<( [ ds j dxdiPe{r5)-\E + i X B\x{n<\i\<2n}dxdn ||(/.||loo || V(/)||l» . 



Observing that, because of (|5.15p . 



diPe{n)-X{n<\£,\<2n} 



1 

< — 
1 



dif\i\^X{n<\i,\<2n} 



LMO,i;Ll(K|)) 



;£n, with £n 0, 



while of course we have for some Cg > 0, 

d(.f3eiTs)-X{n<\(\<2n} 



< an'. 



Therefore, we deduce from the Holder inequality that we have for all e > 0, 

dCl3eirs)^X{n<\^\<2n} ^ 0, 

in L^{0, t; -LP(R^)) for ah 1 < p < | as n ^> cxd, and hence in particular in L^{0, t; L^(R'^)). 
This implies 



1 



f3eiT5)-X{n<\^\<2n}dmdx ^ 0, in ((0, i) X M^) (5.16) 



ds 



as n oo, since E e L°°(0, t;L^( 

On the other hand, using (j2.17p and the fact /3e(r5) < /, we obtain 

1 



ds / dxdC/3e{Ts)-\^\\B\x{n<\^\<2n} 



n 



< - f ds( I mdi] \B\dx 
n Jo \Jr3 J 



< - fds( [ m'dA'\B\dx 

n Jo \Jr3 ) 

< —t sup j [ /ICpdQ ll^llL°°(0,t;L5(R3)) 



(5.17) 



^0, 

as n — )■ oo. Hence, combining (j5.16|) and (|5.17|) together, we get 

[ ds [ dxd£,PeiTs)-\E + £, X B\x{n<\i\<2n} 

as n — >■ oo. 

Finally, letting first n go to oo and then e go to in (|5.14p . we deduce, by Fatou's 
lemma, 



/ Ts{x,^,t)dxd^ <0, for all t>0, 

JR6 
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which, combined with (j5.12p . imphes that = on x (0, oo) almost everywhere. □ 

In other words, 7<5(/"') weakly converges to 75(/). Since 75 is strictly concave on [0, 00), 
we deduce from classical functional analysis arguments that /" converges in measure to / 
on X (0,r) for all T G (0, 00), see |14j . This convergence implies that 

r^f in LP{^,T-L\W^)), (5.18) 

for all 1 < p < 00 and T £ (0, 00). Indeed, by the equi-integrability of the sequence 
{r}n=i and the integrability of / in L^i[0,T] X M^), we know that for any e > 0, there 
exists R > and 6 > such that 

/ \r - f\dtdxdi < e, (5.19) 



sup 

nenJ{[0,T]xBnxBu) 



and 



sup / \r - f\dtdxdi < e, (5.20) 

neN Jg 

for ah set G C [0,r] x with \G\ < 6. 

On the other hand, on the set [0, T] x Br x Br, since up to a subsequence, /" converges 
to / almost everywhere, by Egorov's theorem, for the given 6 > as above, there exists 
a subset H C [0,r] x Br x Br with \H\ < 6 such that /" converges uniformly to / on 
([0,r] X Brx Br) n i?^. Therefore, using (lOoD . 

\r - f\dtdxdi 

[0,T]xBiixBr 

= [ \r-f\dtdxdc+ f \r-f\dtdxd^ (5.21) 

J([0,T]xBRxBR)nH<= Jh 

<e+ [ \r - f\dtdxdi. 

J([Q,T]xBRxBn)r]H'= 

Notice that the last term in (|5.2ip tends to as n — )■ 00 since the uniform convergence of 
/" to / in H^. Hence, combining (j5.19p . (j5.20p and (|5.2ip . we conclude that 

r^f in Li(0,r;Li(M6)), 

which, with the uniform bound of /" in L°°(0, T; L1(IR6)), implies ([5T8]l . 



5.3. Step Three: The convergence in C([0, T]; L^(M^)). It only remains to show that 
converges to / in C([0, T]; ^^(M*^)) using (l5J8]l . Indeed, because of (IXTHI) and g25]), 
it is clearly enough to show that, for each 5 > 0, T € (0, 00), K compact set in M^, we 
have 

/35(r)^/35(/) in C{%T]-L\K)). (5.22) 
For this purpose, we take 4> G C^(M^) such that = 1 on iC, (/> > 0, and we use (j4.3p 
to deduce that for all t > 0, 



i3s{rf<pdxdi= f [ dxdd^^xQirj 



l + '^r (5.23) 
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Then, due to (|538|) . ^siD converges to /3s{f) in LP{M.^ x (0,r)) for all 1 < p < oo and 
T € (0, oo), and one can check easily that the right-hand side of (j5.23p converges uniformly 
in t € [0,T] to the same expression with replaced by /. Since f3s{f) is a renormalized 
solution, this expression is also given by jj^g f3s{f)'^(j)dxd^. In other words, we have 

/ f3s{n^cPdxd(^ [ f3s{ffcpdxd(, (5.24) 

uniformly in t € [0, T], for all T € (0, oo). 
In addition, since ()4.3p implies 

^^^gP e L^(o,r; w-"'1(m6)) 

for large enough n > 0, and 

(3s{n€LH0,T;LHR^)), 

by the Aubin-Lions lemma, we know that /3s{f"^) converges to (3s{f) in C([0, T]; Wj^^'^ {M.^)) 
for any s > 1. Therefore, if we consider L| = L^(supp(^, (pdx), since {/35(/")}n is bounded 
in L°°(0, T;L|), we deduce that (^sif"") converges uniformly on [0, T] to (3s{f) in 
endowed with the weak topology, which, combined with (I5.24p and the fact that f3s{f) € 
C([0, oo); L"^) implies that (^sif^) converges to l^sif) in L| strongly and uniformly in [0, T]. 
Hence, ([02]) fohows. 



6. Large Time Behavior 

In this section, we are devoted to the study of the large time behavior of the renormalized 
solution to VMB. Indeed, let /(t,x,^) be a renormalized solution to VMB with finite 
energy and finite entropy in view of ()3.16p . Then, for every sequence {tnjj^i going to 
infinity, there exists a subsequence {infej^Li ^'^'^ ^ local time-dependent Maxwellian m 
such that fnk{t-,x^£,) = f{t + tn,^,x,^) converges weakly in L-^((0,T) x M^) to m for every 
r > 0. More precisely, we have the following theorem: 

Theorem 6.1. Let /(t,x,^) be a renormalized solution to VMB and assume that b > al- 
most everywhere. Then, for every sequences tn going to infinity, there exists a subsequence 
trif, and a local time- dependent Maxwellian m{t,x,(^) such that fnki'tiX,^,) = f{t + tn,.,x,£,) 
converges weakly in L^{{0, T) xM^) to m{t, x, ^) for every T > 0. Moreover, the Maxwellian 
satisfies the Vlasov-Maxwell equations: 

om 

— + e- V^m + (E + e X B) • V.m = 0, (6.1a) 
at 

dE f 

— - Vx5 = - / m^dS,, divB = 0, (6.1b) 

dB f 

+ V X £; = 0, divE = / mdi, (6.1c) 
dt 7k3 



in the sense of renormalizations. 
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Remark 6.1. When the spatial domain is a periodic box or a bounded domain with the 
reverse reflexion boundary or the specular reflexion boundary, we can expect, as in [71 112]. 
that the local Maxwellian m in Theorem 16.11 is actually global; that is, m is independent 
of t, X. 

Remark 6.2. Our large time behavior result is only sequential; that is, the Maxwellian 
could depend on our choice of the sequence {tnjj^i- 

Proof of Theorem \6.1\ Notice that since is a renormalized solution to VMB, it 

automatically holds: 



sup ( / f{l + \^\^ + v{x) + \\ogf\)dxdi+ [ {\E\^ + \B\^)da 

tG[0,oo) VjR6 JR3 



+ dx didi. dojbiff:-ff,)iogLJ±<^. 

Jo Jr3 Js^ J J* 



(6.2) 



Therefore, fn{t,x,^) = f{t + tn,x,^) is weakly compact in L^((0, T) x M^) for every 
r > and each sequence of positive numbers {tn}^=i going to oo. Similarly, En{t,x) = 
E{t + tn,x), Bn{t,x) = B{t + tn,x) are weakly compact in L°°{0,T; L^{M.^)). Then, the 
weak compactness of fn{t,x,^) in L^((0,T) x M^) implies that there exists a subsequence 
{tnt,}kLi ^-iid a function m E L^{{0,T) x M^) such that the function fn^. converges weakly 
to m in L^[{0,T) x M^) while the weak compactness of Bn{t,x) and En{t,x) implies that 
we can choose tn^. such that and En,, converge weakly* to B and E respectively in 
L°°{0,T;L'^{R^)). Notice that, applying the velocity average lemma, we know 

/ fnj^^ f mdi in L\^,T-L\R?)), 

and 

JR3 Jr3 

Hence, according to (jl.lbp and (jl.lcp . the electric field E and the magnetic field B satisfies 

dE 



V X 5 = - / mCd^, 
with 

div5 = 0, div^;= / md^, 



in the sense of distributions. 

In order to prove that m is a Maxwellian, we denote 

,T r r r f 'f ' 



dk-= dx didi^ I dujb{fn,^'fnJ^-fnJn,,Jlog 

Jo Jr3 Jr6 J52 
T+tnfe r r r ff 



kJ'n-k* 



dx / dS,d(^. / du}h{f f^- f f^)\o, 



J J* 

Then, the estimate (16. 2p implies that dk converges to as A; goes to 00. 
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On the other hand, in view of the first statement of Theorem l3.1l or arguing as [11] , for all 
smooth nonnegative functions ^p, (j) with compact support, we have, up to a subsequence, 



and 



d^di^ / dujbm{t,x,(,')m{t,x,^l)(j){(,)ij{(,^), 



d^d^^ / dub/nJukMOi^i^*) 

d^d^* / dujbm{t,x,^)m{t,x,^^,)(t){C)ij{i*), 
Js2 



(6.3) 



(6.4) 



for almost all {t, x) £ [0, T] x 

Furthermore, since C(M^) is separable, we can also assume the convergence in (j6.3 



and (j6.4p holds for all nonnegative function in C(]R'^). Since P{x,y) = {x — y)ln(|) is a 
nonnegative convex function for x, y > 0, we have, 

0< / d^d^* / dLub{m'm'^ — mm^) log -ip{C*)4'{C) 

Jr6 Js2 mm,:. 



< liminf dk = 0, 

k—^oo 



for almost all {t,x) G [0,r] x R^. Hence, 



b{m'm'^ — mm^) log -^{i')4'{i) = 0) 

mm* 



almost all (t, x) G [0,T] x R'^. The nonnegativity of the function P{x,y) and the strict 
positivity of b ensure that 



m!m! = mrritf^ 



for almost ah (t, x, ^, ^=,,,0;) G (0, T) x x 5^. According to Lemma 2.2 of [3], ] or Section 
3.2 of [5], m is a Maxwellian. Thus, 

Q{m, m) = 0. 

Also, in view of Theorem 13.11 m is still a renormalized solution to VMB, hence 

dim 

— + ^ V^.m +{E + ixB)- V^m = 0, 
at 

in the sense of renormalizations. The proof is complete. □ 
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7. Remark on The Relativistic Vlasov-Maxwell-Boltmann Equations 

An extension of our analysis is possible to the relativistic Vlasov-Maxwell-Boltzmann 
equations of the form (cf . ^ [T7] ) : 

^ + i-'^xf + (.E + ixB)-V^f = Q{f,f), xGM^ CgM', t>0, (7.1a) 
dE 

V X B = -j, divB = 0, on x (0, oo), (7.1b) 



dt 
dB 

'dt 



V X E = 0, divE = p, on x (0, oo), (7.1c) 
P= I fdC, j= [ fidC, on m3^(0,oo), (7.1d) 



with 



and 



Qifj)=[ dcJ du A f'fi-fh). 

Jr3 Js2 ^/l + £ 2^/1 + £ 12 



The corresponding conservation laws are given by 

dp 



+ diva^j = 0, 



d_ 
di 



f^d^ + ExB 



and 



div^ ( / e ® Udi + ( \^^L±\^id -E®E-B®B\]=0, 



d_ 
'di 



y^jVTTWd^+\E{t,x)\^ + \B{t,x)\^'^ 

+ div^ (^^^ m^fd^ + 2E{t, x) X B{t, x)^ = 0. 

Then, we can deduce that for any t G [0, T] 

I f{t,x,i){,/TTW+V^+W)dxdi<C{T), (7.2) 

which implies /|.^| € L°°{(),T;L^{M^). Hence, following the lines in Section 2, the exis- 
tence of renormalized solution to the relativistic Vlasov equation can be verified. More 
importantly, we can further release the requirement on the integrability of the electric field 
in L^, since we no longer need the estimate on in L°°(0,T; L^(]R^)). 

Note that for the relativistic VMB, the magnetic field has the same integrability in the 
variable x as the magnetic field due to equivalence between ^ and y^l + |£p when £ is 
sufficiently large. More precisely, we have 
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Proposition 7.1. Assume that / G L°°((0,T) x M^). Then for any solution satisfying the 
above conservation laws, one has 



where the constant C depends on the energy of the initial data and on ||/||i,°o((o,t)xR'3)- 
Proof. Indeed, we have 

p{t,x)= [ f{t,x,i)di+ I f{t,x,OdC 
<C^\\f\\L^+R-' f fVT+Wd^ 

where for the last inequality, we optimize R by taking 

R = fVl+WdC^ ' . 

The same computation also works for j. □ 

For any sequence of /" as in Section 3, by the H Theorem and (|7.2p . /" is weakly 
compact in L^{{0,T) x (M^). And then, we can follow the lines in Section 4 and Section 5 
to show the corresponding weak stability for the relativistic VMB. One difference is that, 
due to Proposition 17. 11 we need to assume the electric field E(t, x), and the magnetic field 
B(t, x) are uniformly bound in L°°(0, T; L"(M^)) for some a > 4. When the weak stability 
and the existence of renormalized solutions to (|7.ip are concerned, a different assumption 
on the collision kernel need to assume, that is, 

(l + jzp)"M / -^^2=dA^0, as |z|^oo, for all i?G (0,oo). 
\Jz+Bn VI + / 
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